Quantum Blockade and Loop Currents in Graphene with Topological Defects 
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We investigate the effect of topological defects on the transport properties of a narrow ballistic 
ribbon of graphene with zigzag edges. Our results show that the longitudinal conductance vanishes 
at several discrete Fermi energies where the system develops loop orbital electric currents with 
certain chirality. The chirality depends on the direction of the applied bias voltage and the sign of 
the local curvature created by the topological defects. This novel quantum localization phenomenon 
provides a new way to generate a magnetic moment by an external electric field, which can prove 
useful in nanotronics. 
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PACS numbers: 72.10.Fk, 72.15.Rn, 73.20.At, 73.20.Fz, 73.61.Wp 

Graphene is a single layer of graphite with a hon- 
eycomb lattice consisting of two triangular sub-lattices. 
This peculiar structure of graphene gives rise to two lin- 
ear "Dirac-like" energy dispersion spectra around two de- 
generate and inequivalent points K and K' at the cor- 
ner of the Brillouin zone [H, The valley index that 
distinguishes the two Dirac points is a good quantum 
number, even in the presence of weak disorder, since 
inter- valley scattering requires the exchange of large mo- 
mentum. Valley-dependent physics has been actively ex- 
plored recently and can potentially play an important 
role in future graphene based devices d,P3, Q. 

To manipulate this extra degree of freedom, it is nec- 
essary to couple the two valleys. In graphene, there is 
a natural way to produce a valley coupling by creat- 
ing topological defects such as pentagons and heptagons 
Q. Such defects cannot be constructed from a perfect 
graphene sheet simply by replacing a hexagon by a pen- 
tagon or a heptagon. Instead, a "cut-and-paste" process 
0,11] should be employed to keep the local coordination 
number of each carbon atom, as illustrated in Fig. [1] 
A pentagon (heptagon) will induce a positive (negative) 
curvature around it. As can be observed from Fig. [TJ 
after going around any closed carbon loop encircling the 
defect, which always consists of odd number of atoms, 
the roles of two sublattices are interchanged. Therefore, 
the defect breaks the bipartite nature of the lattice in 
the real space as well as the symmetry between the K 
and K' points in reciprocal space, leading to a Mobius 
strip-like structure coupling the two valleys. Theoret- 
ically, the effect of topological defects on the low en- 
ergy electric physics of graphene is equivalent to generat- 
ing non-Abelian gauge potentials with the internal gauge 
group involving the transformation of valley index. The 
wave function acquires a topological phase when circling 
around the defect, which can be described by means of a 




FIG. 1: (color online) The process of constructing a pentagon 
((a)— >(b) — »(c)) and a heptagon ((a)— ►(d)— »(e)) by "cut" and 
"paste" from a perfect graphene sheet. Positive (negative) 
curvature is induced by a pentagon (heptagon) in the realis- 
tic three dimensional space. Notice that after the "cut" and 
"paste" , a single pentagon (heptagon) will give rise to a global 
change of the lattice structure, which is explained in the text. 



non-Abelian gauge field. Experimentally, pentagon and 
heptagon topological defects have been found in graphite- 
related materials [y, Q. In graphene, recently observed 



lfj . 1 1 II ] mesoscopic corrugations (ripples) are partially 



attributed to topological defects. 

The equilibrium electronic properties of a two di- 
mensional graphene in the presence of single or many 
topological defects have already received wide attention 
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However, there is still 
much to understand about the transport properties of 
these systems. In this paper, we investigate the elec- 
tronic transport properties of a zigzag edge graphene 
nanoribbon with several topological defects: a pentagon, 
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a heptagon, and a pentagon-heptagon pair at the center. 
We numerically calculate the total conductance G and 
the spatial distribution of local currents z(r). We reveal 
a novel quantum localization phenomenon whereby the 
conductance vanishes at discrete Fermi energies in the 
first quantized plateau. This effect is accompanied by 
the development of circular loop currents with prescribed 
chirality, owing their existence to the non-Abelian gauge 
potentials connecting the valleys. The chirality depends 
on both the direction of the applied bias voltage and the 
sign of curvature created by topological defects, and can 
hence be easily controlled and manipulated. The result 
opens a new possibility to generate magnetic moments 
by an external electric field in graphene. 

The 7r electrons in graphene are modeled by the tight 
binding spinlcss Hamiltonian 



+*y>k-+H.c), 



(1) 



where c] (a) creates (annihilates) an electron on site i, eo 
(set to 0) is the on-site energy and t (^2.7eV) is the hop- 
ping integral between the nearest neighbor carbon atoms 
separated by a (~1.42A). t and a are used as the energy 
unit and the length unit, respectively. The hopping t is 
taken to be constant for all the C-C bonds, even those on 
the defect. Small changes of the hopping integral t due 
to orbital mixing from the curvature 17| are neglected 
since they do not generate large momentum scattering, 
which means that pure topological effects are considered. 

The zero temperature two terminal conductance G of 
the sample and local density of states (LDOS) p m of site 
m at Fermi energy E F are given by [181 ] 
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where G r{a \E F ) = (E F -H G -T, r s {a) (E F )-T^ a) (Ep))- 1 
is the retarded (advanced) Green's function, G a (E F ) = 
(G r (E F )Y, T S(D) (E F ) = i[Z r s{D) (E F )-Z a siD) (E F )}, H G 

is the Hamiltonian of the graphene sample and Eg = 

H GS G r s {a) HsG = H GS (E F ±i V - H 8 )- l H SG is the re- 
tarded (advanced) self-energy due to the semi-infinite 

source, and the self-energy due to the drain E^> is sim- 
ilar. The self-en ergi es can be numerically calculated in 
a recursive way [19 . l20j ]. We attach clean graphene as 
source and drain leads to avoid redundant scattering from 
mismatched interfaces between the sample and the leads. 

The current per unit energy at the Fermi level E F be- 
tween neighboring sites m and n is HI, 21, 22, 23| 



,(Ef) 



4e 

T 

Ae 



4e 



H mn \ip m \\ip n \ sm(6 r , 




FIG. 2: (color online) Conductance of the topologically dis- 
ordered sample with a pentagon G (black solid), averaged 
loop current {in} (blue dashed dot) and induced magnetic 
moment m (red solid) of the pentagon as functions of Fermi 
energy E. The conductance for a perfect graphene sample Go 
(black dashed) with the same dimensions is also plotted for 
reference. 



where ipm is the amplitude of the wave function at site m, 
G n - G R (T s f{n s )+T D f(fi D ))G A is the electron corre- 
lation function and / is the Fermi distribution function, 
9 n is the phase of %p n and H mn = t. The quantum cur- 
rent i contains information about the phase shift of wave 
functions between the neighboring sites. 

The local current I m ^ n is the energy integral of i m —> n : 



Im — >r 



l (E)dE 



/ — Im[G R (E)T s {E)G A (E)} mn dE I zero temperatu 
ho h 

(ns)eVb I linear response limit 



Im[H mn G™ m (E F )], 



where the bias V& is related to the source (drain) chemical 
potential us (a*d) by ns~ M-D = e ^6- Equation (J7J) can be 
derived from equation ([5]) in a straightforward way along 
with the fact that Green's functions are symmetric ma- 
trices in the absence of magnetic field and thermal fluctu- 
ations, thus the contributions below [Ld cancel, making 
the local current a Fermi surface property [l8j ]. 

We constructed a pentagonal defect in the center of 
a zigzag edge graphene sheet with 45 x 20 sites, in the 
presence of the pentagon, as shown in Fig. [21 the con- 
ductance G is suppressed by defect scattering in most re- 
gions of energy compared to the quantized conductance 

241 of the ballistic case. The G(E) curves exhibit dips 
'aaw kinks at some special energy points. The details of 
these points are the main focus of this work. 

("Bet us first concentrate on the dips A and B in Fig. 
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[2] where the conductance is reduced by the amount 
of go = 2e 2 /h, suggesting a single conducting channel 
at this energy has been completely blocked. The lo- 
cal shape of the dips is approximately Lorentzian. An 
anti-resonance with such features has been observed in 
one-dimensional (ID) systems with topological defects 
2H HI 22 1- The location of these anti-resonances is not 



universal, but is determined by the topological configura- 
tion of the lattice and its defects. Unlike the ID case, the 
locations of these anti-resonances do not have a simple 
relation with the energy levels of the defect or the sam- 
ple. The multi-dimensional nature of the transmission 
matrix and the inter-channel scattering make an analyt- 
ical treatment more complicated. However, the physical 
result is the same as that in ID. The phenomenon at 
hand can be called topological quantum localization. 

To gain better insight in the microscopic origins of the 
localization, we plot, in Fig. [3] (a) and (b), the spatial 
distribution of LDOS p, and local currents i far (point 
J in Fig. [2|) away from and near (point H in Fig. [2|) 
the antiresonance respectively. The LDOS does not vary 
much at the two points, but the currents do. In Fig. [3] 
(a), where G is not sharply suppressed, the spatial dis- 
tribution of i(r) is rather uniform, and almost all of the 
bond currents flow from the source (right) to the drain 
(left). On the contrary, in Fig. [3] (b), loop currents flow 
circularly around the defect in the same direction, giving 
rise to a vortex-like pattern of local currents. Interest- 
ingly, these loop currents reverse their chirality when the 
Fermi energy is swept from one side to the other side 
of the antiresonance, as can be seen in Fig. [3] (c) . The 
magnitude of the current on the pentagon bonds may 
be larger than the source-drain current. Moreover, loop 
currents on some of the hexagons in the proximity of the 
defect can also be observed. This is a quantum effect, 
which is forbidden in a classical resistance network [jjj ]. 
due to the lack of an external battery or magnetic field on 
the loops. This pattern exhibits almost perfect five-fold 
rotation symmetry of the pentagon. Since the graphene 
sample used here lacks the rotational symmetry of the 
pentagon, the pattern in Fig. [3] (b) does not originate 
from an artificial symmetric arrangement of the device, 
but is an intrinsic property of the topological defect, and 
can be expected to exist in the realistic system. In the 
case of an infinite graphene sheet, the first conductance 
plateau containing this bounding state encircling the de- 
fect degenerates into the Dirac point. The existence of 
this bounding state at Dirac point is consistent with pre- 
vious predictions based on the the continuum description 
of 2D graphene @,[l,[l3|. 

The vortex-like loop currents suggest that the topo- 
logical defect introduces an effective magnetic field 29(, 
which is consistent with theoretical predictions [3, H, Hf . 
According to ([5]), the loop currents reflect the pecu- 
liar phase correlation induced by the topological defect. 
Given a closed carbon loop, we define the bond current 




FIG. 3: (color online) The spatial distribution of LDOS of 
the sites (blue circles), local current density of the bonds (ar- 
rows) and induced magnetic moment of carbon rings (green 
disks and crosses) near the pentagon. The size of each symbol 
is proportional to the local value of the corresponding quan- 
tity, (a) E = — O.lt, corresponding to point J in Fig. [5] far 
away from the anti-resonance, (b) E = O.L, corresponding 
to point H in Fig. [5] near the anti-resonance, (c) E — 0.2t, 
corresponding to point Q in Fig. [2] on the other side of the 
anti-resonance. Some circular loop currents around the pen- 
tagon are stressed in purple arrows. Note that circular loop 
currents can also be observed on some hexagons near the de- 
fect. 



on the loop as positive if it points clockwise, and nega- 
tive otherwise. Thus the averaged loop current (ALC) 
density («z,)ioop of a closed loop can be defined as 



(3 L/ loop 



TV 



'■max '-mm 



{k,l) Giloop 



(9) 



where N is the number of bonds (also the number of 



4 



sites) on the loop, 6(x] is the Heaviside step function, 
and i ma ,x (*min) is the maximum (minimum) current on 
the loop. Obviously a current on the loop can be called 
"circular" only if i max ■ i min > 0, so that (iz,}ioo P ^ 0. 
The induced magnetic moment density on a closed loop 
is expressed as [2l| 

m loop (.E)= **^i(-B)(rfc x r,)/2, (10) 

(k,l) Eloop 

where is the coordinates of the site k. A relation 
similar to that of / and i in (JJJ and [5]) holds between 
the magnetic moment M and its energy density m. Due 
to equations (J7J) and ©, we simply discuss the densities 
i and m in the following, but refer to them simply as 
current and magnetic moment. The ALC and magnetic 
moment as functions of energy are plotted in Fig. [21 
They change signs when energy crosses an anti-resonance 
(A and B in Fig. \2§- This is an important correlation 
between the magnetic moment direction and the anti- 
resonances. 

We also note that at the edges of quantized conduc- 
tance plateau (edges of sub-bands also), such as at the 
points D, E, and F in Fig. [21 the conductance has a 
kink-like or a dip-like behavior. This enhanced scatter- 
ing was also observed in disordered graphene with non- 
topological impurities [3(J, and can be attributed to ex- 
treme level-broadening (or velocity renormalization) due 
to van Hove singularity at the sub-band edges [3l|. As 
can be seen from Fig. [21 this enhanced scattering is 
also accompanied by a (qusi-) loop current on the defect, 
albeit with different behavior from the purely topolog- 
ical localization previously discussed: the directions of 

(^L/pentagon 

and/or (m) pentagon are singly-polarized near 
the conductance dip or peak, while the magnitude has a 
sharp peak. 

The transport properties of graphene with a single 
heptagon are also analyzed. FigH] and Fig[5] are the re- 
sults for constructing a heptagonal defect in the center 
of a zigzag edge graphene sheet with 30 x 16 sites. The 
quantum localizations, vortex-like loop currents around 
the defect and reversal of the chirality across the anti- 
resonance can also be observed, except the pattern of 
loop currents now has seven-fold rotationally symmetry. 

We now consider the pentagon-heptagon at the center 
of graphene. With almost zero curvature, this is believed 
to be a stable configuration in carbon- related materials. 
In Fig. the conductance G, the ALCs (iz,) P entagon 
and (ji)heptagon are plotted. Within the first conduc- 
tance plateau, the anti-resonance only happens at the 
center and the edges of the plateau, and (ji)pentagon and 
(«i)heptagon always point in opposite directions. On the 
higher conductance plateaus, we observe different energy 
dependence of the two loop currents: non-zero loop cur- 
rent on the pentagon tends to appear for E > while 
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FIG. 4: (color online) Conductance of the topologically dis- 
ordered sample with a heptagon G (black solid), averaged 
loop current (il) (blue dashed dot) and induced magnetic 
moment m (red solid) of the pentagon as functions of Fermi 
energy E. The conductance for a perfect graphene sample Go 
(black dashed) with the same dimensions is also plotted for 
reference. 

it tends to appear for E < on the heptagon. This 
can be seen more clearly when we compare the local cur- 
rents near the anti-resonance at negative (Fig. [7J (a)) 
and positive energy (Fig. [7J(b)), respectively. There are 
also vortex-like loop currents encircling the defect. When 
E < 0, the magnitude of loop current on the heptagon 
is larger than that on the pentagon, and the loop cur- 
rents encircling the whole defect have the same direction 
with the heptagon loop current (Fig. UJ (a)). Viceversa, 
when E > 0, the magnitude of loop current on the pen- 
tagon is larger, and the vortex-like loop currents have 
the same direction with the pentagon loop current (Fig. 
[7J (b)). All these observations lead to a nontrivial conclu- 
sion that, when the pentagon-heptagon pair is present, 
the pentagon tends to trap the electronic motion in the 
positive energy region, while the heptagon is dominant in 
the negative energy region. The LDOS has a large mag- 
nitude on the defect, similar to the case of disordered 
carbon nanotubes [261 ] . 

The quantum localization is realized only in the first 
quantized channel as the conductance never vanishes 
completely in higher conductance plateau (C and P in 
Fig. [2]). In a graphene narrow ribbon with zigzag edges, 
the first quantized level is different from higher levels due 
to the presence of edge states Q. In the first quantized 
level, for a given Fermi level Ep, the currents carried in 
the states of one valley are rectified currents and the di- 
rections are opposite in different valleys. However, in the 
higher quantized levels, the states in each valley carry 
current in both directions. In the presence of the topo- 
logical defects, the non-Abelian gauge potential scatters 




FIG. 5: (color online) The spatial distribution of LDOS of the 
sites (blue circles), local current density of the bonds (arrows) 
and induced magnetic moment of carbon rings (green disks 
and crosses) near the heptagon. The size of each symbol is 
proportional to the local value of the corresponding quantity, 
(a) E — 0.17t, corresponding to point A in Fig. [4] far away 
from the anti-resonance, (b) E = — 0.15t, corresponding to 
point B in Fig. [4] near the anti-resonance, (c) E = — 0.12t, 
corresponding to point C in Fig. [5] on the other side of the 
anti-resonance. Some loop currents around the pentagon are 
stressed in purple arrows. Note that circular loop currents 
can also be observed on some hexagons near the defect. 



the states in one valley to states in the other valley. In 
the case of the first quantized level, the scattering can 
block the transport current and create the loop current, 
since the process reverse the current completely, while 
this is not the case in the high quantized levels. 

It is interesting to compare our results with those of 
other carbon related structures. In Ref. [21|, a single 
Ceo molecule with discrete energy levels is connected to 



FIG. 6: (color online) Conductance G (black solid) of the 
topologically disordered sample with a pentagon-heptagon 
pair, averaged loop current of the pentagon {it,) pen (blue 
dashed dot) and the heptagon {iz,)hep (red solid) as functions 
of Fermi energy E. 



is two-dimensional (2D) with continuous energy bands. 
Therefore, the role of the degenerate resonance level in 
Ref. 



2l| is quite analogous to that of the anti-resonance 



ID leads, where pentagonal loop currents were discov- 
ered. While, in our work, the sample along with the leads 



level here, where the loop currents reverse their direc- 
tions. The perfect symmetry may be responsible for a 
much larger loop current in Ceo (several tens of Isd) 
than in disordered graphene (few times of Isd)- Elec- 
tronic transport through topologically disordered carbon 
nanotubes (CNTs) has been investigated [26]. But the 
topological defect therein should possess certain configu- 
ration to satisfy the periodic boundary condition in the 
transverse direction. Moreover, the circular motions of 
the electron was verified to play important roles in the 
transport, which is lacked in finite graphene ribbon with 
fixed boundary condition. 

Before arriving at final conclusions, it is important to 
remark that, the understanding of bounding states asso- 
ciated with a single topological defect discussed here can 
be a starting point for investigating the effect of mul- 
tiple topological defects. The interferences between the 
bounding states of the defects may give rise to some rich 
phenomena. 

In conclusion, we have performed a numerical calcu- 
lations of conductance and local currents for topologi- 
cally disordered graphene with one pentagon, heptagon 
or pentagon-heptagon pair. A microscopic understand- 
ing of the conductance reduction is obtained. The strong 
scattering in these systems is always accompanied by 
(quasi-) loop currents around the defects. The chirality 
of the loop currents can be controlled by the bias volt- 
age as well as the gate voltage near the discrete Fermi 
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FIG. 7: (color online) The spatial distribution of LDOS of 
the sites (blue circles), local current density of the bonds (ar- 
rows) near the pentagon-heptagon pair, (a) E = — 0.18t, cor- 
responding to point A in Fig. [6] (b) E = 0.2t, corresponding 
to point B in Fig. [6] Some loop currents near the defect are 
plotted in purple arrows. 

energies where quantum localization takes place. The 
magnetic moments generated by the loop currents can 
be measured in experiments, such as scanning tunneling 
microscope(STM) and Kerr effects. In the presence of the 
pentagon- heptagon pair defect, the pentagon (heptagon) 
is more likely to trap electrons with positive (negative) 
energy. 
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